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ABSTRACT 

The  maximum  strength  of  initially  imperfect,  axially  compressed,  circular 
cylindrical  shells  has  been  studied  with  the  use  of  Reissner's  variation¬ 
al  principle,  von  Iffitrmdn-Donnell  shell  theory,  and  a  deformation  theory  of 
plasticity.  The  results  of  the  present  analysis  reflect  families  of  load- 
end  shortening  curves  for  long  circular  cylinders.  For  a  given  material, 
each  curve  relates  not  only  to  an  imperfection  parameter  which  provides  a 
loading  path  into  the  inelastic  range  but  *»lso  to  the  radius -to-wall 
thickness  ratio.  Significant  maximum  strength  reductions  are  obtained, 
relative  to  predictions  based  on  linear-elastic  theory,  for  specific 
materials  in  the  range  of  radius-to-thickness  ratios  of  practical  inter¬ 
est. 
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INTRODUCTION 


For  almost  forty  years,  many  investigators  have  attempted  to  resolve  the 
discrepancy  between  classical  theory  and  experiment  for  thin  circular 
cylindrical  shells  subjected  to  uniform  axial  compression.  Within  the  last 
several  years,  three  comprehensive  survey  papers  on  shell  buckling  have 
been  prepared  by  Hoff  ’  and  Stein  ,  with  the  strongest  emphasis  directed 
to  thin-walled,  circular  cylindrical  shells  in  axial  compression.  Both 
survey  authors  attribute  the  discrepancy  between  the  theory  and  experiment 
to  the  combined  effects,  in  varying  degrees,  of  prebuckling  deformations, 
initial  imperfections,  and  boundary  conditions;  yet,  except  for  noting  the 
possible  limitation  of  linear  elastic  theory,  the  quantitative  effect  of 
inelastic  deformations  on  the  initial  buckling  (maximum)  strength  of  such 
shells  is  neither  discussed  nor  estimated  by  either  of  the  authors  on  the 
basis  of  the  literature  surveyed.  However,  in  view  of  the  experimental 
work  of  a  number  of  investigators  directed  to  initially  imperfect,  axially 
compressed,  circular  cylindrical  shells  (see  refs.  4,5,6,and  7),  not  only 
have  inelastic  deformations  been  detected  at  buckling  (maximum  load)  but 
they  have  been  shown  to  reduce  the  load-carrying  capacity  significantly  as 
a  result.  Further,  several  approximate  analyses  (refs.  4,7,8,  and  9)  have 
indicated  the  presence  of  inelastic  deformations  at  initial  buckling  of 
imperfect  shells.  (The  inelastic  buckling  of  imperfect  shells  investigated 
by  Lee^  is  noted  here  but  not  discussed  further,  as  interest  is  being 
centered  on  shells  which,  in  the  absence  of  inperfections,  would  buckle 
only  in  the  elastic  range.)  Since  the  circular  cylindrical  shell  is  a  basic 
structural  element  in  aerospace  design  and  since  its  initial  buckling  load 
is  indeed  its  maximum  load,  a  clear  understanding  of  all  of  the  factors 
influencing  it6  initial  or  maximum  load  under  axial  compression  and  the 
relative  magnitudes  of  these  factors  is  necessary. 

Now,  in  the  case  of  plates  compressed  beyond  the  classical  buckling  load, 
it  is  well  known  that  a  unique  load-shortening  curve  exists  for  all  plates 
of  a  given  material  independent  of  width-to-thickness  ratio,  is  long  as  the 
material  follows  a  linear-elastic,  stress-strain  law.  However,  a  series  of 
analyses  by  Mayers  et  al.  ’  ’  ,  leading  to  an  essentially  exact  solution 

for  the  maximum  strength  of  compressed  rectangular  plates,  based  upon 
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von  K&rnrin's  large-deflection,  strain-displacement  relations,  has  shown 
that  every  plate  of  a  given  material,  hut  different  width-to-thickness 
ratio,  compressed  beyond  the  classical  buckling  load  into  the  Inelastic 
range,  possesses  a  unique  load-end  shortening  curve  and,  hence,  maximum 
load.  Similarly,  existing  von  IfttrTn£n-Donnell-type  elastic  postbuckling 
analyses  of  axially  compressed,  long,  thin  cylinders  indicate  again  a 
unique  load-shortening  relationship  for  all  radius -to-thickness  ratios. 
Now,  as  a  result  of  the  present  study,  families  of  load- shortening  curves 
have  been  obtained  for  initially  imperfect,  long  circular  cylinders. 

Each  curve  for  a  given  material  relates  to  a  particular  value  of  the 
radius -to-wall  thickness  ratio  and  to  an  initial  imperfection  pattern 
which  provides  .1  loading  path  into  the  inelastic  range,  since  the  classi¬ 
cal  buckling  load  for  compressed  cylinders,  unlike  plates,  is  never 
reached  in  reality.  In  the  range  of  radius-to-thickness  ratios  of 
practical  interest  for  both  unstiffened  and  stiffened  shells,  significant 
maximum  strength  reductions  are  obtained  relative  to  predictions  based  on 
any  theory  utilizing  linear  stress-strain  behavior. 

The  maximum  strength  behavior  of  initially  imperfect,  axially  compressed, 

circular  cylindrical  shells  has  been  studied  with  the  use  of  a  modified 

14 

form  of  Reissner's  variational  principle  and  the  von  K&rm&i-Donnell 
strain  displacement  relations.  In  general,  Reissner's  variational  prin¬ 
ciple  not  only  permits  the  selection  of  the  stresses  independent  of  the 
displacements  but  also  facilitates  the  incorporation  of  inelastic  effects 
into  the  analysis.  Confidence  in  the  inelastic  analysis  based  on 
Reissner's  principle  has  been  established  herein  by  comparing  the  special 
case  of  the  purely  elastic  solution  with  that  obtained  by  Kempner1'’ 
through  the  use  of  the  minimum  total  potential  energy  principle.  A  sim¬ 
ilar  procedure  was  followed  with  respect  to  plates  in  the  analyses  of 
reference  13.  With  confidence  thus  established  in  the  method  of  solution, 
inelastic,  maximum  strength,  radius-to-thickness-dependent,  load-shorten¬ 
ing  curves  have  been  obtained  for  shells  of  different  materials.  To 
obtain  the  present  maximum  strength  criteria  for  thin  shells,  the  problem 
i A" investigated  qualitatively  and  economically  by  using  a  beam-arch  analog 
for  the  cylinder  and  quantitatively  by  analyzing  the  circular  cylindrical 
shell  on  the  basis  of  analog-generated  trends. 
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GENERAL  THECM 

STATEMENT  OF  THE  PROBLEM  AND  BASIC  ASSUMPTIONS 


The  problem  considered  is  that  of  the  maximum  strength  analysis  of  a 
long,  initially  imperfect, circular  cylindrical  shell.  The  solution  for 
the  shell  maximum  strength  is  obtained  through  the  use  of  Reissner's 
variational  principle  modified  to  include  the  effects  of  nonlinear- 
elastic  material  behavior.  No  unloading  in  the  nonlinear-elastic  range 
is  assumed  to  occur;  thus,  the  terms  nonlinear-elastic  and  Inelastic  are 
used  interchangeably.  The  presence  of  unloading  in  the  plaBtic  range  is 
investigated  a  posteriori, and  the  effect  of  unloading,  if  any,  on  the 
solution  obtained  is  assessed  in  the  discussion  of  the  results  of  the 
analysis.  To  establish  confidence  in  the  inelastic  analysis,  the  elastic 
solution  for  the  postbuckling  problem  is  extracted  and  compared  with  that 
developed  on  the  basis  of  the  minimum  total  potential  energy 
principle . 

As  In  the  plate  analyses  of  references  11-13,  a  two-element  description 
of  the  shell  cross  section  (Figure  l)  is  used  to  avoid  the  complexity,  due 
to  inelastic  effects,  of  integrating  a  nonlinear  stress  distribution 
through  the  cylinder  wall.  The  core  separating  the  two  elements  is 
assumed  to  be  rigid  in  shear,  and  the  inelastic  behavior  (constitutive 
equations  and  yield  criterion)  is  based  upon  simple  deformation  theory. 
Material  compressibility  effects  in  the  inelastic  analysis  are  assumed 
to  be  negligible;  therefore,  Poisson's  ratio  Is  taken  as  one-half  when¬ 
ever  the  load-shortening  relationship  reflects  nonlinear  behavior. 

The  von  K&rraAn-Donnell  strain-displacement  equations  for  an  initially 
imperfect,  thin-walled,  circular  cylindrical  shell  are  used  in  the 
present  analysis.  For  the  radial  deflection-thickness  ratio  magnitudes 
involved  in  an  initial  buckling  or  maximum  strength  analysis  of  a  thin 

o 

shell  in  compression,  and  in  view  of  the  findings  of  Mayers  and  Rehfield 
based  on  a  more  accurate  nonlinear  theory  for  thin  shells,  the  use  of  the 
von  K^rmAn-Donnell  strain-displacement  relations  is  deemed  justified. 
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BASIC  EQUATIONS 

The  von  Kdrmifn-Donnell  strain-displacement  relatione  for  the  circular 
cylindrical  shell  are  modified  to  include  the  effect  of  initial  imperfec¬ 
tions.  Initial  imperfections,  as  defined  herein,  are  deviations  from  the 
exact  circular  cylindrical  shape  of  the  shell  which  are  present  before 
any  axial  compression  takes  place. 

The  total  displacements  of  the  shell  can  be  expressed  as 
Uj  *  u 


W_  =  V  +  W 

T  o 

where  w  is  the  initial  radial  deviation  of  the  shell  from  its  exact 
o 

shape  and  u,  v,  and  w  (see  Figure  l)  are  the  additional  displacements  of 
the  shell  from  the  initial  shape  during  the  loading  process.  The  mid- 
surface  strain-displacement  equations  in  terms  of  total  displacements  are 
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VARIATIONAL  PRINCIPLE 

The  Reissner  variational  principle  for  prescribed  surface  displacement 
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(shell  under  controlled  end-shortening)  selects  from  among  1  states  of 
stress  and  displacement  which  satisfy  the  boundary  conditions  of  prescribed 
surface  displacement,  the  true  state  of  stress  and  displacement  by  requir¬ 
ing  the  satisfaction  of  the  equation 

6U  ’  6{jffir<Vx+V y  +  WF',dV}  =0  (5) 


where  the  stress  and  displacement  states  are  varied  independently.  F', 
the  stress  energy  density,  is  a  function  of  the  stresses  such  that  the 
stresses  and  strains  are  related  by  the  constitutive  equations 
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For  linear  problems,  F'  becomes  the  complementary  energy  density.  The 
principle  establishes  the  Euler  equations,  which  consist  of  three  equi¬ 
librium  equations  and  three  stress-displacement  relations,  as  well  as  the 

associated  boundary  integrals  from  which  natural  boundary  conditions, 
if  any,  can  be  established.  The  Euler  equations  and  boundary  conditions 
for  the  circular  cylindrical  shell  behaving  according  to  the  von  Kdrmrin- 
Donnell  theory  are  developed  in  Appendix  I. 


The  stress  energy  density  F'  is  defined,  in  general,  as 
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From  the  secant -modulus  deformation  theory  of  plasticity,  for  an  incom¬ 
pressible  material,  the  "effective"  strain  can  be  written  as 
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and  "effective"  stress  as 
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where 


eff 


~  Eseeff 


and  E  is  the  secant  modulus  of  the  uniaxial  stress-strain  curve  of  the 
s 

material.  For  the  present  analysis,  the  uniaxial-stress  strain  curves  in 
tension  and  compression  are  taken  as  identical.  Based  upon  equation  (j), 
the  incremental  form  of  F'  can  be  written  as 
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As  shown  in  reference  11,  the  incremental  quantity  dF'  can  be  expressed 
also  in  the  form 
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Therefore, 
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The  Reissner  functional  is  then  expressed  as 
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With  the  introduction  of  the  Ramberg-Osgood  three -parameter  representation 

of  a  uniaxial  stress -strain  curve  in  the  form 
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the  functional  U"  becomes 
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Consistent  with  the  two-element  description  of  the  shell  and  after  inte¬ 
gration  over  the  thickness  h,  the  result,  as  developed  in  Appendix  I,  is 
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where  tu  =  x/L,  e  =  y/2«R,  and  V  =  4jtt^.LR. 
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The  total  stresses  and  strains  are  written  in  terms  of  midsurface  and  bend¬ 
ing  contributions, as  denoted  by  the  primed  and  double-primed  quantities. 

The  subscript  t  refers  to  the  top  face  and  the  subscript  b  refers  to 
the  bottom  face  of  the  two-element  cross  section.  With  the  von  Kdrman- 
Donnell  strain-displacement  equations  and  the  relations  for  the  effective 
stress  in  terms  of  the  primed  and  double-primed  quantities,  as  expressed 
in  reference  13,  the  Reissner  functional  becomes 
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In  consideration  of  the  linear-elastic  problem  for  a  homogeneous  cylinder, 
and  with  the  bending  stresses  related  to  the  curvatures  through  Hooke's 
law,  a  modified  Peissner  function  (following  ref.  13)  is  expressed  in  the 
form 


This  modified  form  of  Reissner's  principle  leaves  only  the  displacements 
and  midsurface  stresses  as  the  variational  quantities.  As  described  in 
the  next  section,  equation  (l8)  is  used  as  the  basis  of  correlation  with 
Kempner's  minimum  potential  energy  solution  (ref.  15)  for  the  postbuckling 
load-shortening  curve. 


METHOD  OF  SOLUTION 

Previous  investigators  (Cox"^,  Tsien^®,  and  Hoff^)  have  approximated  the 
postbuckling  behavior  of  both  axially  compressed  and  uniformly  heated, 
circular  cylindrical  shells  by  analyzing  various  one -dimensional  models. 
These  models  consist  of  a  basic  beam  member  supported  by  one  or  more 
springs  whose  load-deflection  curves  are  chosen  such  as  to  force  the  models' 
load-shcrtening  curves  to  provide  the  snap-through  behavior  associated  with 
the  pos  ibuckling  behavior  of  thin  shells  in  compression. 


In  the  present  analysis,  a  model  consisting  of  a  beam  resting  on  a  series 
of  shallow  sinusoidally  curved  arches,  as  shown  in  Figure  2,  is  selected 
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to  be  the  analog  of  the  shell.  The  behavior  of  an  arch  under  radial 
loading,  shown  in  Figure  3,  provides  a  nonlinear  foundation  for  the  beam 
so  as  to  yield  the  snap-through  behavior  associated  with  axially  com¬ 
pressed,  cylindrical  shells.  The  beam  length  is  equivalent  to  the  cylinder 
length;  the  arch  length  is  taken  as  a  circumferential  wave  length  of  one 
buckle.  The  spacing  between  arches  is  equivalent  to  the  axial  wave  length 
of  one  buckle  of  the  cylinder.  The  radius  R  of  the  cylinder  is  brought 
into  the  analog  geometry  through  the  rise  parameter  for  circulaily  curved 
plates  of  small  chord;  that  is,  the  maximum  rise  aQ  of  the  sinusoidally 
curved  arch  in  Figure  3  is  equated  to  the  rise  of  a  shallow,  circularly 
curved  plate  through  the  relationship  aQ  =  A^/8r . 

SHALLOW  ARCH  BEHAVIOR 


A  Reissner  functional  for  a  shallow  arch  subjected  to  a  prescribed  uniform 
loading  is 


U"  +  V"  r'AC  (S'  ,  0  u  or" 

"Hi - /  {/  [v’y  +  2V’y~p]  +  E~[2  W’yy] 

0  V 


where  V"  m  potential  of  the  prescribed  loading 

—  =  variable  curvature  =  -  w 

P  o,yy 

wq  =  initial  shape 

p  =  uniform  load 

DISPLACEMENT  AND  STRESS  FORMULATIONS 

Functions  for  the  displacements  and  stresses  are  chosen  in  the  form 
w  a 

r  =  r  6ln  ^ 


(19) 
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The  use  of  a  vanishing  v-displacement  in  conjunction  with  the  modified- 
Reissner-principle  approach  and  the  von  KdrmAn  strain-  displacement 

equations  has  been  justified  in  reference  13  for  plates.  For  the  arch, 
with  the  kinematics  described  by  a  specialization  of  the  von  Karman- 
Donnell  curved  plate  equations,  the  stresses  and  w  displacement  can  be 
determined  from  the  first  variation  of  the  Rei6sner  functional  regardless 
of  the  Fourier -series  representation  for  v.  Thus,  in  the  absence  of  any 
possible  coupling  of  cr  and  v  in  the  strain  energy,  it  is  obvious  that 

y 

v  can  be  equated  to  zero  without  any  loss  in  generality. 

LOAD-DEFLECTION  CURVE 

Substitution  of  equation  (20)  into  equation  (19),  with  subsequent  integra¬ 
tion  over  the  arch  length  and  variation  with  respect  to  the  amplitudes  of 
the  displacements  and  stresses,  results  in  three  simultaneous  equations  for 
the  determination  of  the  nondimensionalized  displacement  and  stress  param¬ 
eters  for  a,  o^,and  cr^  defined  in  equation  (20).  The  simultaneous 
equations  can  be  written  conveniently  in  terms  of  only  one  of  the  param¬ 
eters  (equation  (60)  of  Appendix  II);  that  is, 


3  3£ 

1  "  2t| 


,2  t  («  +  l)t  .  <i£)2  1 


where 


I  =  a/h 


11  ~  2  a 


5  -  S<t> 

A  plot  of  p  vs.  |  is  shown  in  Figure  3. 
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It  is  interesting  to  note  that  the  single  variationally  derived  equation 

above,  for  the  case  of  a  linear-elastic  constitutive  law  (Hooke's  law), 

reproduces  exactly  the  load-deflection  curve  for  a  shallow  arch  constructed 
20 

by  Timoshenko  from  two  equations  derived  from  beam  theory  with  midplane 
stretching  effects  included. 

BEAM-ARCH  ANALOG 

The  Reissner  functional  for  the  entire  beam-arch  model,  including  the 
effects  of  inelastic  behavior  and  the  use  of  the  two-element  ’ascription 
for  the  beam  and  arch  cross  sections,  is  written  as 


UA  UB  ffjcf  1  2  1 
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where  <p  =  y /Aa,  u)  =  x/^. 

DISPLACEMENT  ARP  STRESS  FORMULATIONS 

In  addition  to  equation  (20),  expressions  for  the  beam  displacement  and 

stresses  are 
u  =  -ex 
O' 

x  o 
E  =  "  E 


(23) 


where  e  =  prescribed  end  shortening.  The  use  of  midsurface  displacements 
involving  no  free  parameters  in  conjunction  with  the  modified-Reissner- 
principle  approach  has  been  justified  for  plates  and  shells  in  references 
9  and  13  and  is  demonstrated  herein  for  the  arch.  Thus,  the  stresses  and 
w  displacement  can  "be  determined  from  the  first  variation  of  the  Reissner 
functional  regardless  of  any  additional  Fourier-series  representation  for 
u;  therefore,  it  is  sufficient  to  describe  u  as  -ex. 

LQAD-SHQRTEMIHG  RELATIONSHIP  (ELASTIC  PROBLEM) 

For  elastic  behavior,  K  is  set  equal  to  zero  in  equation  (22).  As  shown 
in  Appendix  II,  substitution  of  equations  (20)  and  (23)  for  the  assumed 
stresses  and  displacements  of  the  arch  and  beam,  respectively,  into  the 
total  Reissner  functional  (eq.  (22))  for  the  analog,  followed  by  integra¬ 
tion  over  one  wave  length  in  each  direction  and  variation  with  respect  to 
«n  free  parameters,  leads  to  five  simultaneous  equations  which  can  be 
reduced  by  elimination  to  only  two  equations;  that  is, 


These  equations  are  solved  in  the  following  manner.  First,  for  given  values 
of  the  parameters  eR/h,  p,  r|,  and  Bq,  equation  (25)  is  evaluated 
numerically  for  £.  Then,  for  the  given  parameter  values  and  the  corre¬ 
sponding  i's,  equation  (24)  is  solved  numerically  to  establish  the 
relationship  between  <JR/Eh  and  eR/h;  namely,  the  family  of  imperfection- 
dependent  load-shortening  curves  shown  in  Figure  4. 
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LOAD-SHORTENING  RELATIONSHIP  (INELASTIC  PROBLEM) 

The  inelastic  behavior  is  studied  by  retaining  the  two  terms  multiplied 
by  the  parameter  K  in  equation  ( 22 ).  To  enter  the  inelastic  range,  an 
elastic  analysis  is  repeated  at  an  end -shortening  value  which  insures 
elastic  behavior.  The  parameter  values  obtained  are  utilized  as  starting 
points  for  the  inelastic  analysis.  A  Newton-Raphson  procedure  is  employed 
to  establish  the  load-shortening  curves  for  inelastic  behavior  with 
increasing  values  of  the  end  shortening.  Families  of  radius-to-thickness- 
dependent,  load-shortening  curves  are  then  obtained  for  E,  K,  and  n 
values  corresponding  to  electroformed  nickel  (see  stress-strain  curve  in 
Figure  5)-  Such  load-shortening  curves  are  presented  in  Figure  6.  The 
thickness  t  of  the  homogeneous  material  is  related  to  the  distance 
separating  the  two  faces  of  the  model  by  setting  the  critical  load 
of  the  two-element  shell  equal  to  the  critical  load  of  the  homogeneous 
shell  of  given  R/t.  The  resulting  relationship  between  h  and  t  is 
simply  h/t  =  l/\fz. 

CIRCULAR  CYLINDRICAL  SHELL 
DISPLACEMENT  AND  STRESS  FORMULATIONS 

The  maximum  strength  of  a  long, circular  cylindrical  shell  is  studied  for 
both  the  elastic  and  inelastic  problems  by  enforcing  the  vanishing  of  the 
first  variation  of  functionals  given  by  equations  (l8)  and  (17),  respec¬ 
tively.  The  variations  are  carried  out  with  respect  to  assumed  states  of 
stress  and  displacement,  with  the  latter  being  kinematically  admissible. 

The  approach  followed  is  that  of  prescribing  the  stresses  and  the  radial 
deflection  pattern  in  the  same  form  as  those  formulated  by  Kempner'*'^. 
After  a  point  on  Kempnsr's  elastic  postbuckling  curve  is  reproduced  by 
the  elastic  results  of  the  present  analysis  (K  =  0  in  eq.  (22)),  confi¬ 
dence  in  the  elastic  analysis  is  established.  The  radial  displacement 
in  any  one  buckle  of  the  shell  is  selected  to  be 
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The  corresponding  midsurface  displacements  would  in  general  be  taken  as 


OD  00 


However,  the  use  of  expressions  for  u  and  v  involving  no  free 
parameters  in  conjunction  with  the  modified-Reissner-principle  approach 
has  been  demonstrated  in  references  9  and  13  and  discussed  earlier  with 
reference  to  both  the  arch  and  the  beam-arch  analog.  It  has  been  shown 
that  the  stresses  and  radial  displacement  may  be  determined  from  the  first 
variation  of  the  Reissner  functional,  regardless  of  the  values  of  u^ 
and  v  when  the  von  Kdrman -Donne 11  strain-displacement  relations  are 

1  J 

used  along  with  Fourier-series  representations  for  both  displacements  and 

stresses.  Thus,  It  is  once  more  concluded  that  u  and  v  can  be 

1 J  1 J 

equated  to  zero  without  any  loss  in  generality;  the  midsurface  displace¬ 
ments  become  simply 

u  =  -ex  (27) 

v  =  0  (28) 

The  rationale  underlying  the  utilization  of  midsurface  displacements 
reflecting  no  free  parameters  in  references  9  and  13  as  well  as  the  present 
work  is  simply  the  realization  that  either  nonexistent  or  very  weak  cou¬ 
pling  exists  between  free  midsurface  stresses  and  displacements  in  a 
Reissner  formulation, depending  upon  the  degree  of  nonlinearity  in  the  plate 
or  shell  kinematics  being  employed. 

The  simplification  of  the  expressions  for  the  midsurface  displacements 
Involving  no  free  parameters  permits  the  equilibrium  equations  in  the  x 
and  y  directions  to  be  satisfied  independently  of  the  magnitudes  of  the 
free  stress  coefficients. 

The  inplane  stresses,  in  addition  to  the  average  axial  stress,  consist  of 
the  Fourier  harmonics  found  in  Kempner's  stress  function  (ref.  15)  and  are 
expressed  in  the  form 

ll* 
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It  is  noted  that  the  variable  coefficients  in  a '/E,  t'  /e,  and  o'/E 

x  3ty  y 

are  the  same.  This  equivalence  is  established  automatically  from  the 
satisfaction  of  the  inplane  equilibrium  equations  in  the  x-  and  y-directions, 
respectively,  prior  to  employing  the  variational  principle,  because  the  dis¬ 
placements  u  and  v  in  equations  (27)  and  (28)  are  not  arbitrary. 

Now,  instead  of  relating  the  bending  stresses  to  the  curvatures  through 
Hooke's  law,  they  are  selected  in  the  form  of  the  curvatures,  but  with 
arbitrary  coefficients,  and  are  written  as 
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(32) 


(33) 

(34) 


The  Justification  for  this  assumption  i6  established  in  references  9  and 
13  through  comparison  of  results  with  those  of  well-known,  essentially 
exact  solutions  for  postbuckling  analyses  of  shells  and  plates,  respec¬ 
tively. 
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Because  of  the  periodic  nature  of  the  assumed  expressions  for  the  nonuniform 
stresses  and  displacements,  the  boundary  integrals  derived  on  the  basis  of 
the  Reissner  variational  principle  in  Appendix  I  can  be  reduced  to 


for  any  y 
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(37) 


It  can  be  seen  immediately  that  expressions  (29),  (30),  and  (3l)  for  the 
stresses  satisfy  the  integral  conditions  (35),  (36),  and  (37). 

INITIAL  IMPERFECTION  PATTERN 
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In  a  recent  report  by  Tennyson  and  Welles  ,  the  authors  note  that  the 

first  terms  to  be  observed  experimentally  at  the  inception  of  buckling  are 

those  described  by  the  combination  of  an  asymmetric  plus  an  axisymmetric 

term.  Moreover,  from  their  analysis  of  free  vibrations  of  thin-walled 

22 

cylinders,  Mayers  and  Wrenn  point  out  that  for  finite  displacements  of 
the  order  of  magnitude  of  the  shell  thickness,  the  minimum  energy  is 
related  to  a  slightly  "modified-chessboard"  radial  deflection  pattern 
given  by  the  same  harmonics  as  those  in  reference  21.  It  is  conjectured, 
therefore,  that  in  view  of  the  random  nature  of  imperfections  in  thin 
shells,  a  reasonable  imperfection  shape  to  which  the  shell  might  be  partic¬ 
ularly  sensitive  is 

wo  eh  [_hiQ 

LOAD-SHORTENING  RELATIONSHIP  (ELASTIC  PROBLEM) 
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To  establish  the  all-elastic,  load-shortening  curves  for  the  circular 
cylindrical  shell,  a  stationary  value  of  the  Reissner  functional  U"  with 
respect  to  the  free  parameters  g  ,  A^,  and  0/ E  is  sought  (see  Appendix 
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Ill)*  The  set  of  equations  established  from  the  vanishing  of  the  first 
variation  is  reducible  to  four  equations  in  the  four  unknowns  cf/E, 

and  Igg*  This  set  of  four  equations  is  then  used  to  construct  the 
imperfection-dependent,  elastic  load-shortening  curves  presented  in 
Figure  7  Tor  n  =  1  and  q  =  0.45*  These  equations  are  used  also  to 
reproduce  a  point  on  the  Kempner  load-shortening  curve  for  ji  =  0.3899* 
q  =  0.226,  eR/t  =  0.4,  and  &=  0.0  as  a  means  of  establishing  the 
accuracy  of  the  analysis. 

LOAD-SHORTENING  RELATIONSHIP  (INELASTIC  PROBLEM) 

In  the  inelastic  problem,  the  load-shortening  curve  for  the  circular 
cylindrical  shell  is  obtained  by  finding  the  extrema  of  the  following 
functional: 
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which  consists  initially  of  17  free  parameters. 


The  additional  free  parameters  over  the  four  required  for  the  elastic 
solution  appear  because  the  bending  stresses,  even  though  taken  in  the 
form  of  the  curvatures,  require  seven  additional  free  parameters  when  the 
deflection  function  is  as  given  by  equation  (26),'  further,  the  six  mid- 


surface  stress  components,  A 


i  y 


cannot  be  eliminated  in  terms  of  the 


i^j  and  o/E  because  the  inelastic  contribution  to  the  total  Reissner 
functional,  which  contains  the  and  a/ E,  requires  that  the  extrema 

of  equation  (39)  be  found  numerically.  As  a  result,  equation  (39)  is 
written  in  terms  of  the  17  free  parameters  51,52,53,A11,A22,A;l3,A13,A02, 

A20 ' &11  * a20 * a02 *^11 ' b20 ’ b02 ; ^11  ’  and  a/E*  7316  double  integral  in  equa¬ 
tion  (39)  is  evaluated  numerically  by  a  Newton-Raphson  iteration  technique, 
and  maximum-strength,  radius-to-thickness-dependent,  load -shortening 
curves  are  obtained  for  three  kinds  of  materials:  2024-T3  aluminum, 
electrof ormed  nickel,  and  electroformed  copper.  Except  for  the  radius- 
to-thickness  ratio,  the  parameter  values  used  are  the  same  as  those  taken 
for  the  elastic  analysis  which  is  applicable  to  all  values  of  R/t.  A 


unique  load-shortening  curve  (independent  of  R/t)  for  all-elastic 
analyses  is  typical  of  solutions  based  on  the  von  K^rm^n  theory  for 
plates  and  the  von  K^rm^n-Donnell  theory  for  shells  (see,  for  example, 
refs.  1,  2,  3,  and  13). 

RESULTS  AND  DISCUSSION 

The  maximum  strength  of  initially  imperfect,  axially  compressed,  circular 
cylindrical  shells  has  been  studied  through  the  use  of  modified  forms  of 
Reissner's  variational  principle,  in  conjunction  with  the  von  Karmdn- 
Donnell  theory  and  a  deformation  theory  of  plasticity.  Maximum- strength, 
radius-to-thickness-dependent,  load-shortening  curves  have  been  obtained 
for  three  different  materials, utilizing  both  a  beam-arch  analog  to  represent 
the  cylinder  and  a  two-dimensional  analysis  of  the  circular  cylindrical 
shell  itself.  The  fundamental  materials  considered  are  nickel  and  copper, 
in  view  of  the  fact  that  most  of  the  experimental  work  carried  out  in 
recent  years  to  study  imperfection  sensitivity  relative  to  axial  compression 
loading  of  carefully  manufactured  (electroformed)  shells  has  dealt  with 
these  materials.  Results  are  presented  also  for  2024-T3  aluminum  shells, 
since  design  practice  reflects  such  a  more  conventional  material. 

BEAM-ARCH  ANALOG 

With  the  use  of  the  beam-arch  analog,  maximum -strength,  radius-to-thickness- 
dependent,  load-shortening  curves  have  been  obtained  that  represent  cylin¬ 
ders  with  material  properties  corresponding  to  electroformed  nickel.  For  a 

given  imperfection  amplitude,  the  load-shortening  curves  for  electroformed 
nickel  reflect  a  radius-to-thickness  dependence  of  significant  magnitude,  as 
depicted  graphically  in  Figure  6.  The  upper  curve  represents  the  results 
of  an  all-elastic  analysis.  The  family  of  curves  belov;  the  purely  elastic- 
analysis  curve  is  obtained  by  incorporating  inelastic  effects  into  the 
variational  principle,  using  the  stress-strain  curve  of  Figure  5  and  per¬ 
turbing  from  the  original  elastic-solution  parameters  established  in 
developing  the  all-elastic  behavior  curve.  It  is  noted  that  for  R/t  =  300 
and  A  =  0.1,  the  reduction  in  maximum  load  from  the  maximum  load  obtained 
in  the  purely  elastic  analysis  is  1 6$.  Also,  for  this  arbitrary  case,  it 
may  be  seen  that  for  a  sufficiently  "thin"  shell,  say,  R/t  >  1000,  the 
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inelastic  analysis  reflects  the  absence  of  inelastic  strains  and  coalesces 
with  the  unique  elastic  analysis  result  (all  R/t). 

CIRCULAR  CYLINDRICAL  SHELL 

For  circular  cylindrical  shells,  maximum- strength,  radius-to-thickness 
dependent,  load-shortening  curves  have  been  obtained  utilizing  stress- 
strain  curves  representing  several  different  materials.  The  survey  papers 
of  Hoff  *  and  Stein  indicate  that  a  number  of  investigators  have  achieved 
unusually  high  buckling  stresses  relative  to  the  classical  stress  by 
performing  exacting  experiments  on  very  carefully  fabricated  specimens. 

The  results  of  these  investigators  are  included  in  Figure  8,  along  with  a 
band  encompassing  the  myriad  tests  performed  prior  to  i960.  The  source  of 
the  band  is  an  indication  of  the  significant  discrepancy  and  scatter  that 
existed  between  experiment  and  classical  theory  over  a  period  of  about 
30  years. 

Above  the  upper  bound  of  the  scatter  band,  results  are  shown  for  the  very 

23  5 

specialized  tests  carried  out  by  Babcock  and  Sechler  ,  Horton  and  Durham  , 
21  24 

Tennyson  and  Welles  ,  and  Almroth  et  al.  The  Babcock  and  Sechler  tests 
were  performed  on  electroformed  copper  cylinders,  the  Horton  et  al.  and 
Almroth  et  al.  tests  were  performed  on  electroformed  nickel  cylinders, 
and  the  Tennyson  and  Welles  tests  were  performed  on  a  series  of  spun-cast 
plastic  cylinders.  In  each  instance,  relatively  small  cylinders  were 
involved,  and  the  control  of  any  irregularities  in  the  shell  geometries  is 
obvious  due  to  the  fact  that  the  classical  load  is  approached  within  10 
percent  on  one  end  (R/t  =  100)  and  25  percent  on  the  other  (R/t  =  1000), 
which  is  well  above  the  upper  bound  of  the  scatter  band  for  much  larger 
shells  fabricated  in  accordance  with  common  shop  practices  over  several  decades. 

From  Figure  8,  it  would  appear  that  the  Almroth  et  al.  and  Babcock  and 
Sechler  electroformed  nickel  and  copper  shells  are  the  most  imperfection 
sensitive;  thus,  these  particular  materials  are  of  analysis  interest.  The 
uniaxial  stress-strain  curve  utilized  in  the  analysis  of  the  electroformed 
nickel  shells  has  been  presented  previously  as  Figure  5,  in  conjunction 
with  the  inelastic  beam-arch  analog  leading  to  the  load-shortening  curves 
shown  in  Figure  6.  Now,  with  this  stress-strain  curve,  the  cylinder 
analysis  yields  the  load -shortening  relationships  presented  in  Figure  9» 
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The  results  show  that  for  p  ■  1  and  A  =  0.1,  the  reduction  from  all¬ 
elastic  behavior  to  the  R/t  *  300  cylinder  Is  about  25  percent.  It  is 
noted  then  that  (l)  inelastic  deformations  can  be  present  in  cylinders 
which  if  perfect  would  buckle  at  the  classical  elastic  load  and  that  (2) 
the  phenomenon  of  inelastic  deformations  produces  a  significant  R/t-depen- 
dence  of  the  load-shortening  curve  for  a  cylinder  of  a  given  material  and 
imperfection  parameter. 

The  family  of  load-shortening  curves  shown  in  Figure  9  is  related  to  only 
one  value  of  this  imperfection  parameter.  However,  for  usual  allowable 
design  tolerances,  imperfections  for  cylinders  corresponding  to  large 
radius-to-thickness  ratios  would  be  expected  to  be  "effectively"  greater 
than  those  of  cylinders  corresponding  to  relatively  small  radius-to- 
thickness  ratios.  Therefore,  the  radius-to-thickness  ratios  are  corre¬ 
lated  to  reasonable  imperfection  parameters  (namely,  the  average  values 
of  the  scatter  band  for  practically  fabricated  shells  of  Figure  8)  to 
establish  the  imperfection  parameter  as  a  function  of  R/t;  the  ensuing 
load- shortening  curves  are  shown  in  Figure  10. 

Each  set  of  curves  for  various  R/t  values  in  Figure  10  is  related  to  a 
given  imperfection  parameter;  the  imperfection  parameter  increases  as  the 
radius-to-thickness  ratio  increases.  Clearly,  as  should  be  expected, 
cylinders  corresponding  to  large  radius-to-thickness  ratios  (for  example, 
R/t  ■  800)  possess  lower  maximum  loads  than  cylinders  of  significantly 
smaller  radius-to-thickness  ratios  (for  example,  R/t  =  400) .  This 
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phenomenon  has  been  demonstrated  experimentally  by  Kanemitsu  and  Nojima 
and  Weingarten  et  al.^  it  may  be  seen  from  Figures  9  and  10  that  as  the 
radius-to-thickness  ratio  decreases,  the  effect  of  inelastic  deformations 
increases.  In  Figure  10,  the  load-shortening  curves  for  R/t  =  400,  600, 
and  800  depart  significantly  from  the  elastic  curves  for  the  same  imper¬ 
fection  amplitude,  ranging  from  about  8  percent  for  the  thinnest  cylinder 
to  25  percent  for  the  thickest  cylinder,  which  is  still  relatively  thin 
by  design  standards.  It  may  be  noted  also  that  cylinders  of  different 
geometries  may  exhibit  the  same  maximum  load  as  shown  in  Figure  10.  A 
cylinder  of  R/t  ■  600  and  A  =  0.1  has  essentially  the  same  maximum  load  as 
a  very  thin  cylinder  (that  is,  R/t  >  1000)  and  A  =  0.16.  The  maximum  load 
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of  the  first  cylinder  depends  not  only  on  the  initial  imperfection  ampli¬ 
tude  but  also  on  local  yielding  of  the  material,  whereas  the  maximum 
load  of  the  second  cylinder  depends  only  upon  the  relatively  large  imper¬ 
fection.  This  is  attributed  to  the  fact  that  bending  predominates  and 
prevents  the  development  of  membrane  strains  which  are  of  sufficient 
magnitude  to  combine  with  the  bending  strains  and  cause  local  yielding 
of  the  material.  As  a  result,  the  total  effective  stress  remains  well 
below  the  yield  point  of  the  uniaxial  stress-strain  curve  for  the  materi¬ 
al,  and  essentially  elastic  behavior  governs. 

Now,  it  is  essential  to  study  more  closely  the  actual  nickel  and  copper 
cylinder  tests  reported  by  Almroth  et  al.  and  Babcock  and  Sechler  for 
cylinders  of  about  R/t  =  850  and  890,  respectively.  The  stress-strain 
curve  for  electroformed  copper  is  shown  in  Figure  11;  the  corresponding 
stress-strain  curve  for  electroformed  nickel  appears  in  Figure  5»  With 
the  fit  of  the  appropriate  imperfection  parameters,  the  maximum  strengths 
produced  by  the  cylinder  inelastic  analysis  are  shown  in  Figure  12.  The 
results  are  obtained  on  the  basis  of  equal  buckle  wavelengths  (p  =  1) 
and  imperfection  amplitudes  (t  =  0.02).  When  it  is  assumed  that  the 
materials  behave  in  linear  elastic  fashion,  the  cylinder  analysis  predicts 
maximum  strengths  some  10-14  percent  higher.  In  other  words f  the  rela¬ 
tively  high  ratio  of  o/oci  reached  in  the  careful  tests  on  carefully 
fabricated  cylinders  could  have  been  higher  and  closer  to  the  classical 
theory  results  reported  by  Stein  ,  had  the  materials  behaved  in  perfectly 
elastic  fashion. 

It  is  believed  that  sufficient  theoretical  evidence  has  been  given, 
in  consideration  of  the  results  of  the  recent  experimental  references 
presented,  to  reach  the  conclusion  that  inelastic  deformation  should 
not  be  excluded  from  initial  buckling  or  maximum  strength  analyses. 

However,  electroformed  nickel  and  copper  are  not  typical  of  practical 
construction  materials.  Thus,  attention  is  next  directed  to  aluminum, 
the  most  common  material  utilized  in  aerospace  construction.  For 
analysis  purpose,  the  stress-strain  curve  of  2024 -T3  aluminum  alloy 
(Figure  13)  is  utilized.  The  resulting  load- shortening  curves  are 
shown  on  Figure  14.  Surprisingly,  the  behavior  of  2024 -T3  aluminum 


cylinders  in  the  maximum  strength  analysis  is  markedly  different  from  that 
of  the  electroformed  copper  and  nickel  cylinders.  The  separation  of 
the  load-shortening  curves  as  a  function  of  R/t  does  not  begin  until 
the  R/t  ratio  reaches  the  neighborhood  of  200.  In  fact,  appreciable 
reduction  in  maximum  strength  does  not  occur  until  a  radius-to-thickness 
ratio  of  about  100  is  approached.  Obviously,  there  is  an  effect  of  the 
nature  of  the  stress-strain  curve  on  the  sensitivity  to  inelastic  effects 
at  a  given  R/t  ratio.  Aluminum  and  other  conventional  materials  exhibit 
much  sharper  knees  of  the  stress-strain  curve,  and  this  can  be  noted 
quickly  from  the  exponent  in  the  Ramberg-Osgood  stress-strain  curve 
representation.  Whereas  the  electroformed  copper  and  nickel  stress- 
strain  curves  reflect  exponents  of  3,  aluminum  and  most  other  conventional 
aerospace  metals  reflect  exponents  of  9  or  greater.  One  exception  is 
stainless  steel;  the  uniaxial  stress-strain  curve  is  represented  quite 
accurately  by  an  exponent  of  3. 

Little  comfort  can  be  taken  from  the  fact  that  the  2024-T3  aluminum 
shells  exhibit  purely  elastic  behavior  in  the  range  R/t  >  300.  Regard¬ 
less  of  choice  of  conventional  structural  material,  thin  shells  for 
resisting  compression  loads  are  utilized  either  in  orthotropically 
stiffened  or  sandwich  forms.  As  pointed  out  by  Hoff  ,  such  shells 
possess  "effective"  R/t  ratios  as  small  as  50.  Hence,  maximum 
strength  determination  for  such  shells  should  be  based  on  analyses 
such  as  the  one  described  here,  which  can  account  for  the  presence  of 
inelastic  deformations.  It  should  be  noted  that  for  "effectively" 
thick  shells,  initial  imperfections  would  tend  to  be  minimized  and  the 
effects  of  inelastic  deformation  maximized.  In  reported  compression 
tests  on  large  design-scale  sandwich  and  orthotropically  stiffened  shells 
(refs.  26,  27,  28),  in  the  range  (R/t)effective>  45-220,  maximum  calcu¬ 
lated  stresses  indicate  inelastic  deformations  at  maximum  load.  The 
family  of  curves  presented  in  Figure  14  illustrates  the  reduction  in 
maximum  load  due  to  plasticity  of  "effectively"  thick  compressed  shells 
of  2024-T3  aluminum;  the  uniaxial  stress-strain  curve  is  shown  in 
Figure  13.  For  R/t  =  60,  the  reduction  in  maximum  strength  due  to 
inelastic  considerations  is  about  33$.  It  should  be  further  noted  that 
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the  imperfection  parameter  value  (A  =  O.l)  would  probably  be  smaller 
in  practice  for  shells  of  effective  R/t  a  100;  thus,  it  would  be  expected 
that  at  higher  average  stress  level  ratios,  o/acl,  ^e  percent  of  reduc¬ 
tions  in  maximum  strengths  due  to  inelastic  effects  should  be  greater  than 
those  shown  in  Figure  14.  A  particular  curve  of  interest  in  Figure  14 
is  that  for  R/t  =  l6j.  This  load- shortening  curve  represents  the  analyt¬ 
ical  prediction  for  a  2024-T3  aluminum  test  cylinder  of  NASA  (reference 

28)  with  the  cross  section  shown  in  Figure  14.  The  imperfection 
parameter,  A  =  0.1,  is  taken  from  the  average  value  of  the  scatter  band  of 
Figure  8  for  R/t  =  167.  In  this  case,  little  effect  of  plasticity  is 
seen.  However,  for  R/t  =  167  and  a  2024-T3  stiffened  aluminum  cylinder, 
the  buckling  mode  is  still  representative  of  thin-cylinder  behavior. 

A  photograph  of  the  buckled  test  cylinder  (see  Figure  15)  clearly  shows 
the  "diamond"  buckle  pattern  typical  of  thin  cylinder  behavior  in  the 
presence  of  initial  imperfections.  An  efficient  stiffened  2024-T3 
aluminum  cylinder  would  appear  to  be  one  of  lower  R/t  than  167,  with 
resulting  smaller  imperfection  amplitude  and  Increased  sensitivity  to 
inelastic  deformations. 

Finally,  for  the  results  obtained  herein,  based  upon  a  simple  deformation 
theory  of  plasticity  with  no  unloading  mechanism  provided,  the  occurrence 
of  unloading,  if  any,  within  a  buckle  must  be  assessed.  In  the  case  of 
each  shell  analyzed,  monitoring  the  numerical  calculations  for  the 
effective  strain  at  a  point,  as  the  end  shortening  increased  incrementally, 
indicated  that  no  detectable  unloading  in  inelastic  regions  of  a  buckle 
occurred  as  the  maximum  load  of  the  shell  was  approached.  Therefore, 
the  initial  assumption  that  the  present  maximum  strength  analysis  could 
be  undertaken  on  the  basis  that  the  material  is  of  a  nonlinear  elastic 
nature  was  indeed  valid. 

CONCLUDING  REMARKS  ' 

The  study  of  the  maximum  strength  of  initially  imperfect,  axially  com¬ 
pressed,  thin-walled,  circular  cylindrical  shells  has  been  undertaken  to 
determine  the  effects  of  inelastic  deformations  on  maximum  load-carrying 
capacity  to  establish,  as  was  found  to  be  the  case  of  compressed  flat 
plates,  the  degree  of  geometry  dependence  of  load-shortening  curves  for  a 
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given  material.  Indeed,  the  load- shortening  curves  obtained  in  the 
present  analysis  for  the  inelastic  behavior  of  shells  reflect  signifi¬ 
cant  radius-to-thickness  ratio  dependence.  A  family  of  load-shortening 
curves,  each  corresponding  to  a  different  value  of  radius-to-thickness 
ratio,  is  shown  to  result  for  a  given  material  and  a  specific  initial 
imperfection  amplitude.  The  significant  spread  between  the  load-shortening 
curves  in  any  family  of  such  curves  is  a  new  result  in  application  of  the 
von  Karman-Donnell  strain-displacement  relations,  which  sheds  further 
light  on  a  problem  whose  satisfactory  solution  has  eluded  investigators 
for  decades.  The  actual  solution  has  been  carried  out  with  the  use  of 
the  von  Karman-Donnell  formulation  and  a  deformation  theory  of  plasticity 
in  conjunction  with  a  modified  form  of  the  Reissner  variational  principle. 

In  view  of  the  results  obtained  for  shells  corresponding  to  experimentally 
studied  electroformed  copper  and  nickel  shells  of  relatively  high  geomet¬ 
ric  accuracy,  it  must  be  concluded  that  initial  imperfections  cannot 
be  considered  the  only  significant  cause  of  such  shells’  failing  to 
develop  their  classical  buckling  stresses.  Further,  the  presence  of 
inelastic  deformations  and  the  nature  of  the  nonlinearity  of  the 
material  stress-strain  curve  for  a  given  shell  cause  considerable  scat¬ 
ter  in  maximum  strength  results  for  cylinders  of  the  same  R/t  ratio. 

Thus,  it  is  reasonable  to  conclude  that  inelastic  deformations,  along 
with  inital  imperfections,  prebuckling  deformations, and  boundary  condi¬ 
tions  (factors  previously  established  to  be  important),  must  be  considered 
in  determirir?  the  maximum  strength  of  thin  shells,  especially  those 
with  small  t  moderate  initial  imperfections.  Such  imperfections  would 
be  representative  of  orthotropically  stiffened  and  sandwich  shells  in 
aerospace  vehicle  designs.  The  more  efficient  stiffened  or  sandwich 
shell  would  appear  to  be  one  in  which  elastic  buckling  would  take  place 
except  for  the  presence  of  some  degree  of  imperfection  which  precipitates 
bending  with  inelastic  deformations.  These  two  effects,  coupled  with 
prebuckling  deformations  and  boundary  conditions,  should  make  optimization 
procedures  more  meaningful. 

A  beam-arch  analog  developed  to  represent  a  circular  cylindrical  shell 
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and  the  shell  itself  have  been  used  to  establish  the  load  shortening 
curves  presented.  The  analog  is  introduced  to  establish  qualitative 
trends  and  to  achieve  economy  in  computation  for  what  must  be  considered 
to  be  an  extremely  complicated  problem  in  the  nonlinear  analysis  of  thin 
shells . 
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Figure  2.  Beam-Arch  Analog  for  Circular  Cylindrical  Shell.  (Pure 
Diamond  Nodal  Pattern  Shown  for  Illustrative  Purposes 
Only.  Shell  Analysis  Is  Based  on  Ovalized  Deflection 
Pattern . ) 
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Figure  3.  Load  Deflection  Relationship  for  Snap-Through  Buckling  of 
Shallow  Arch. 
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Figure  5 •  Stress-Strain  Curve  for  Electroformed  Nickel  Utilized  in 
Present  Inelastic  Analyses. 
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Figure  6.  Family  of  Radius-to-Thickness-Dependent  Load-Shortening 

Curves  Corresponding  to  a  Particular  Imperfection  Parameter 
for  Electroformed  Nickel  Beam-Arch  Analogs. 
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Figure  7 
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.  Load-Shortening  Curves  for  a  Circular  Cylindrical  Shell 
Obtained  for  Elastic  Behavior. 
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Figure  9.  Family  of  Load-Shortening  Curves  for  Electrof ormed  Nickel 
Circular  Cylindrical  Shells. 
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CIRCULAR  CYLINDRICAL  SHELL 


Figure  10.  Families  of  Radius -to -Thickness -Dependent,  Load-Shortening 
Curves  Corresponding  to  Different  Imperfection  Parameters 
for  Electroformed  Nickel  Circular  Cylindiicsl  Shells 
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Figure  11.  Stress-Strain  Curve  for  Electroformed  Copper  Utilized  in 
Present  Inelastic  Analyses. 
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Figure  15-  Buckled  Stringer  Stiffened,  Circular 
Cylindrical  Shell. 
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APPENDIX  I 

EULER  EQUATIONS  AND  BOUNDARY  CONDITIONS  IGRIVED  FROM  THE 
REISSNER  FUNCTIONAL  FOR  A  CYLINDER  WITH  PRESCRIBED  END  SHORTENING 

The  Reissner  functional  in  the  absence  of  prescribed  surface  loading  is 
defined  as 


U"  =  (a  e  +  a  e  +  t  y  -  F'  )dV 
JJJ  x  x  y  y  xyrxy 

V 


The  strain-displacement  relations  of  von  Karman -Donnell.  theory,  modified 
for  the  two-element  cylinder  of  Figure  1  and  including  initial  radial 
deformation  effects  are 


e  =  e 1  ±  e"  =  u,  +  i  w,  +  w  w,  ±  ^  w, 
xt^b  x  x  x  2  x  o,x  'x  2  XX 


€  =  € 1  db  €  =  V  W,  +  W  W,  +  —  w 

yt,b  y  y  ,y  R  2  'y  o,y  'y*  2  'yy 


=  Y'  ±  y"  =  u,  +  v,  +  w,  w,  +  w  w,  +  w  w, 
xyt,b  ^  xy  y  x  X  y  o,x  ;y  o,y  ;x 


±  hw, 


The  stresses  in  terms  of  midsurface  and  bending  components  are 


a  =  o'  ±  or" 


a  =  a'  ±  a" 

yt,b  y  y 


T  =  T*  ±  t" 

xyt;b  *y  *y 


For  the  two-element  section,  equation  (4o)  becomes 
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L  2*R 


u'  =  tf//  fw 

0  0  l  1  t 


ff  e  +ffe  +  o  e  +t  y 
W  yt  yt  yb  yb  xyt  xyt 


+  Txybrxyb  -  (Ft  +  Fb}  )  ^ 


Then,  in  view  of  equations  (4l)  and  (42),  equation  (43)  may  be  written  as 


L  2jiR 


U  =  t 


+  °xbj  [_U*X  +  2  W*X  +  W°>xW,xJ 

+  (%  +  \)  [  ^  +  5  “'y  +  "o,y”-y  ‘  5  ] 

+  /t  +  1  \  r  u»  +  v,  +  W,  W, 

[  xyt  xybj  ’y  'x  'x  'y 

+  w  w  ,  +w  w,  1  +  /  a  -  cr  \  T  ^  W,xx  1 

o,x  y  °,y  x  J  (  xt  xb)[2  J 

i\  -  \  )[  5  “’»]+  (X  ■  X)  [ta'*y]  -(FtX) 


The  average  stresses  and  the  bending  moments  are 


o'  =  ^(o  +  a  ) 
x  2  xt  *b 


o'  =  i(o  +  0  ) 

y  2  yt  yb 

t  1  -  — ( <r  +  t  ) 

xy  2k Txyt  xyb 

M  =  -t.  |(o  -  cr  )  =  -t  ho" 
x  f  2  x,  x,  f  x 

t  b 

M  =  -t_  $(o  -  o  )  =  -t  ho" 
y  f  2  yt  yb  f  y 


M  =  t.  ^(t  -  t  )  =  t  hi" 

xy  f  2V  xyt  xyb  f  xy 
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Substitution  of  equation  (45)  into  equation  (44)  gives 

12.  _  v 


M  M 

.  JLWj  + 

2t  w,xx  2t  ’yy 


x  'y  o,x  'y  o,y  x 

dxdy 


M  ^ 

r^w,  -F'J 
tf  xy  J 


(46) 


With  cognizance  taken  of  the  fact  that 
6F’  =  60  +  6o  +  ^ —  &T 


or 


8F'  =  &cr'  +  8a'  +  5t 

+  8M  +  8M  +  6M 

aiT  x  aiT  y  aiT^  xy 

then  the  vanishing  of  the  first  variation  of  u"  with  respect  to 
V  °y»  ^xy’  \>  My>  Mxy’  Uj  v>  and  w  re<*ulres  that 


6U"  =  2t, 


L  2nE 

/  /  (  8o'(u,  +  ~  wf  +  w,  w  ) 

JJ  \  xK  ’x  2  ’x  ’x  o,x' 
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_  ,  /  ,  1  2  J  w  \ 
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y  y  2  ’y  y  o,y  R 
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xy'  ’x  y  'x  V  o,x  ’y  o,y  x 

+  a'(8u,  +  w,  6v,  +  w  6w,  ) 

x  x  x  x  o,x  x 

\  5w 

+  a'(6v,  +  w,  8w,  +w  8w,  )-a'  — 
yv  y  *y  y  o,y  y  y  R 
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xyv  ’x  ’y  ’x  y  ’x  y 


4  V 


JttUOtk 
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._£w,  - 

-2, 

2tf  xx 

2tf 

M 

M 
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5M 
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uf 
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xy  ^ 


The  result  of  carrying  out  the  integration  by  parts  is 

T.  PitR 


6U"  =  2t. 


(  .  1  2  . 

(u,  +  *  w,  +  w,  w 


dF' 


x'  2  "■'x'  ”,x"o,x  do^B°x 
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y>y  xy,x  x  x  o,x  ,x 

+(c'[w,  +  w  ])  +(t*  [w,  +  w  ]) 

y  y  °>y  >y  xy  ’y  o,y  ,x 

a' 
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Moment -curvature  relations: 
1 


dF' 

ctfT 


dF' 

w 

y 

dF’ 


xy 


2tf  W,xx 


2t„  W,yy 


—  V, 

tf  xy 


(50) 


Midsurface  eQuilibrium  eauations: 


o'  +  t' 
x,x  xy,y 


=  0 


+  t' 


=  0 


y,y  xy,x 
Lateral  equilibrium  equation: 

(°;[w,x+  v0iXl),x  +  (ay(w»y+  wo,y^,y 

+<TJy["’y+  ”o,y]),x  +  <Tiy[",x+  V]l,y 

O’  , 

+  JL  +  _i_fM  +  M  -  2M  3=0 
R  2tfl  x,xx  y,yy  xy,xyJ 

In  view  of  equation  (5l),  the  lateral  equilibrium  equation  (52) 
becomes 

a'(w,  +  w  )  +  a'(w,  +  w  ) 

x'  xx  o,xx'  y'  yy  o,yy' 

a' 

+  2t*  (w,  w  )  + 
xyv  xy  o,xy  R 

+  sr-(M  +  M  -  2M  )  =  0 
2tf  x,xx  y,yy  xy,xy' 


(51) 


(52) 


(53) 


which,for  a  homogeneous,  isotropic  body  (2t^,-»t)  and  Hookean  material, 
reduces  to  the  well-known  lateral  equilibrium  equation  of  von  Karman- 
Donnell  theory  (see,  for  example,  ref.  15). 

The  midsurface  stress-displacement  boundary  integrals  are 


50 
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_ 

f  o'6u| 

'I  X  0 


dy  =  0 


(54a) 
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f  T1  6v| 

J  *y  o 


dy  =  0 


f  a'8v| 

J  y  i 

rL  2nR 
/t1  Bui 

o 


2«R 

I  dx  =  0 


d^  =  0 


The  remaining  boundary  integrals  are 
2*R 

2t„  /  {  a'(w,  +  w  )  +  r'  (w,  +  w  ) 

fj  \  x'  f x  o,x  xyv  y  o,y 

+  M  -  t—  M  )  6w  |  dy  =  0 
2tf  x,x  tf  xy,yj  ^  ^ 


t_  /  ((Ww,  +  w  )  +  t'  (w,  +  w  ) 
f.J  ^y'  y  o,y  xyv  'x  o,x' 


1  i  ^  2J 

-i-  M  -  =-  M  >  6w| 
2tf  y,y  tf  xy,x  J  £ 


dx  =  0 


L  2flR 

2M  Bw|  |  =  0 

^  00 


(54b) 


(54c) 


(54d) 


(55a) 


(55b) 


(55c) 


For  the  present  problem,  that  of  a  long,  circular  cylindrical 
shell,  the  end  shortening  is  prescribed;  thus  Bu  vanishes  at  x  *  0,L. 
The  remaining  boundary  conditions  at  x  =  0,L  can  be  ignored  due  to  the 
length  assumed  for  the  shell.  The  boundary  conditions  at  y  =  0,  2nR 
are  either  continuity  conditions  (55a,  b,  and  c)  or  conditions  which 
preclude  the  transfer  of  any  circumferential  and  shear  loads  across  a 
generator  (54  c,  d). 
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APPENDIX  II 

SHALLCW-ARCH  AND  BEAM-ARCH  LOAD-SHORTENING  RELATI ONSEIPS 
SHALLOW  ARCH 

The  expressions  for  the  displacements  and  stresses  [equation  (20)]  can  he 
substituted  into  the  Reissner  functional  given  by  equation  (19)*  Upon 
subsequent  integration,  the  Reissner  functional  for  a  shallow  arch  under 
uniform  loading  is  then 


VE(h/r) 

a 


<J.R 

A 


_ _ —  H_(!  _  —  | ) 

2  Eh  5  q 

SiWW*  -  #)2 


1(C1A  R  \  2  - 

'12V  E  hi  '  « 


(56) 


where 


I  =  a/h 

4rt2Rh  «2/h 


R  = 


— (— ) 

2  'a„' 


5  ■  gtf’2 

The  first  variation  of  with  respect  to  0 a and  5  results  in 

<r.R  _  /  ^  _2 


gl ar 
Eh 


and 


gAR  Tj/  £\  gLAH  n  2  - 
"  Eh  lS^21  "  q  J  "  Eh  4H'jtP=° 


These  simultaneous  equations  can  be  written  in  terms  of  |  alone  as 
-2  „  I  h  A  /  -wC\2  — 


I 


-S- '’-(if  •«'■(«* 


(57) 

(58) 

(59) 

(60) 
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For  prescribed  values  of  r)  and  p,  equation  (60 )  can  be  solved  numerically 
for  g.  The  results  are  shown  in  a  plot  of  p  vs.  |  in  the  lower  portion 
of  Figure  3. 


BEAM-ARCH  ANALOG 


The  expressions  for  the  displacements  and  stresses  [equations  (20)  and 
(23)]  ca,n  be  substituted  into  the  Reissner  functional  given  by  equation 
(22) .  The  result  of  integrating  the  Reissner  functional  for  a  two-element 
beam  arch  is 


TT,,  PB  1  UA  .  SB/  eR  iVl  SV*  t 

"  ~ ' a(  h  ^  ' 

axR  .  2 
“  Eh" 

0,  ,Bi 


V  1H S  .  1(2?.)"-  &L Ui2~i\ 

2h  ^6  2'Eh'  ^Eli  /  \Eh  JlSV  i\  5/ 

^lA711  R  1  1  /qiAR\2 

"  “eI5”  h  "  2 'Eh  )  “  Jf\  Eh  ' 


(61) 


where 

t  =  a/h 

5  =  a  /h 

o  o' 

y  = 

e  =  applied  unit  end  shortening 

The  first  variation  of  U*  with  respect  to  a,  o  ,  {,  ofl,  and  o  results 

J-  n  In 

in 


(62) 
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hU 
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(63) 
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(continued) 
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tfAR  r\  JZ  AR  « 

sr(i5)(21  ■  r} "  ~  h  =  0 


The  above  equations  can  be  written  in  terms  of  the  parameter  J 
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•l1 -ft -TVS'] 

+  '[-X?  +  ^+ir)  +  1*7] 


-  4r-° 

H  1 


For  prescribed  values  of  n,  q,  6q,  and  eR/h,  equation  (67)  can 
solved  numerically  for  5.  Then,  for  the  same  parameter  values 
corresponding  j's,  equation  (62)  can  be  solved  numerically  to 
the  load- shortening  relationship  between  aR/Eh  and  eR/h  and  to 
the  curves  shown  in  Figure 


(66) 

alone  as 
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APPENDIX  III 

MSTHOD  OF  SOLUTION  FOR  ELASTIC  CIRCULAR  CYLINDRICAL  SHELL 

After  the  expressions  for  the  assumed  displacements  [equations  (26),  (27), 
and  (28)]  and  stresses  [equations  (29),  (30),  and  (31)]  are  substituted 
into  the  Reissner  functional  given  by  equation  (18),  subsequent  integra¬ 
tion  yields 


~  Et(2*RL)  “  E  F0C“  Ao2F02"  A11^FL1+ 

•**221*22+  H22)"3A3l(F31+  3H31> 
"A13^3F13+  H13^‘8A20H20 


-2(f)  -  |k4l+  l64)(l^}  + 

f  a|1(1+9u2)  +  32(A^2+  44a|0)] 

+  ^M2[(1+h2)2|2  +  32(mJ+|2  +  |2)j+4  |  e 
2Ml-v2)  1  2  3  E 


2s2.  ,2 


where  the  F's  and  H's  are  defined  as 


_  2  t 11  .  ,2  ,^11  ,  520n 

Foo =  ^ n  R[~r +  *20 +  2|n  (_r +  ~r)] 

0 


F02  ~  ^  ^  R^ll  2^11  *11^ 

O 


P  f  s02 

F11  =  ^  11  R^2^11^20+^  11^02 +2^ll  ^20 +  ~2~ 

o 

»«- A  «4i  +  *W«n  (1t"¥)! 


_  2  tr2.  ,  /  ^20  .  6lln 

F3i  =  M  ^  r^h^20  2|n  ^3  12  ^ 

o 
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F13  “  ^  71  R^11^02  +  511  ^02^ 

0 


i- 


5- 


u  2  trbll  .  .  5Uv  2  t  . 

H20  =  M  n  r[“5"  +  hi  TT)_M  R  620 

o 


H11  =  4  11  R^lll02+  ^Uq^02  ^  R  ^11 
£2  I 

„  2  tr  511  5ll! 

H22  =  M  11  r  ”  TT  ’  4n  T] 

o 


H31-° 

2  t 

®13  =  M  11  R^ll502-  *11  ^02^ 

o 

In  the  above,  the  substitution  £2Q0  =  In  A  133,8  teen  made  for  conve¬ 
nience  based  on  the  postbuckling  analysis  results  presented  in  reference 

29*  It  is  of  interest  to  note  that  the  same  assumption  has  been 
utilized  by  Tennyson  et  al.  in  reference  21  in  conjunction  with  an 
elastic  analysis  of  the  buckling  process  for  compressed  cylinders. 

The  equations  found  from  the  vanishing  of  the  first  variation  can  be 
reduced  to  four  equations  in  the  four  unknowns  Iq2>  ancl  a/E • 

This  set  of  nonlinear  algebraic  equations  can  now  be  solved  numerically 
for  given  values  of  the  parameters  A,  \i,  and  r|,  and  the  relationship 
between  dR/Et  and  eR/t  can  be  obtained.  Curves  of  crR/Et  vs.  eR/t  are 
given  in  Figure  7» 
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